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We discuss dynamical pairing correlations in the context of configuration mixing of 
projected self-consistent mean-field states, and the origin of a divergence that might 
appear when such calculations are done using an energy functional in the spirit of a 
naive generalized density functional theory. 



1. Introduction 

Self-consistent mean-field models are one of the standard approaches in nuclear 
structure theory, see Refffl for a recent review. For heavy nuclei, they are the only 
microscopic method that can be systematically applied on a large scale. 

Over the last few years, we were involved in the development of a method, that 
adds long-range correlations to self-consistent mean-field by projection after varia- 
tion and variational configuration mixing within the generator coordinate method 
(GCM). A pedagogical introduction to our method is given in we will sketch 

some of our recent work on those aspects of this method that concern the treatment 
of pairing correlations. 



2. Theoretical framework 

2.1. Self- consistent mean-field with pairing 

For an introdu ction i nto the treatment of pairing correlations in a HFB framework, 
see for 

examplJSHS! 

and references given therein. We will give only those details 
that are relevant for the further discussion. 

All results given below were obtained with the effective Skyrme interaction 
and a density-dependent pairing interaction with a soft cutoff at 5 MeV 
above and below the Fermi energy as described in The HFB equations are 
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solved using the two-basis method introduced in Ref A All we need to know about 
an HFB state in what follows is that in its canonical basis it is given by 

|HFB) = nK+«M4«).)|0)- (1) 

in terms of occupation probabilities and creation and annihilation operators. 

HFB states are not eigenstates of the particle-number operator. In condensed 
matter, for which HFB theory was originally designed, this is not much of a problem, 
as the number of particles is usually huge. In nuclear physics, where the particle 
number is quite small, this causes two problems: On the one hand, the standard 
HFB treatment artificially breaks down when the density of single-particle levels 
around the Fermi energy is below a critical value, leading to a HF state without 
pairing correlations. On the other hand, a HFB state mixes the wave functions 
of different nuclei as it is spread in particle-number space, with a width that is 
proportional to the dispersion of the particle numbei!^. 

3. Particle- number projection 

variation after projection provides a rigorous solution to both problems. 

For var ious mainly technical reasons, however, this approach is not yet widely used 
(seJU^for an exception). We use a different approach instead, that was first outlined 
m In a first step, we complement the HFB equation with the Lipkin-Nogami 

(LN) method, that provides a numerically simple approximation to projection be- 

o 

fore variation, see Ref.° and references given therein for details. The LN procedure 
enforces the presence of pairing correlations also in the weak pairing regime at small 
level density. The LN method gives a correction to the total energy whose quality 
has been repeatedly questioned. We do not make use of the correction term, but 
project in a second step on exact eigenstates of the particle-number operator N, 
with eigenvalue A'o, applying the particle- number projection operator 

^^^0 = 7^ rd^e^^(^-^o}^ (2) 



2tt _ ^ 

The operator e'^'^^ rotates the HFB state in a U{1) gauge space 

|HFB(^)) = e^'^^lHFB) - l[{u,, + v^e'^^alal)\0), (3) 

while e""^^" is a weight function. For states with even particle number Nq the 
integration interval can be reduced to [0,7r]. We discretize the integrals over the 
gauge angle with a simple i-point trapezoidal formula 

"^^ ^ 1=1 

As was shown by Fomenkcff^, this simple scheme eliminates exactly all compo- 
nents from an HFB state which differ from the desired particle number A^o by up 



I 



I 
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to ±2(L — 1) particles. Although the spread of the near-Gaussian distribution of 
particle numbers is large in comparison with respect to the constrained value, it is 
small enough that already small values of L, ranging from 5 in light nuclei to 13 in 
heavy ones, are sufficient for the numerical convergence of the integrals over ip. 

In practice, we constrain the HFB states in the mean-field calculations to the 
same (integer) particle number that we project on afterwards. This is, however, 
not a necessary condition. In a projection-before-variation approach, the particle 
number of the intrinsic state has no physical meaning, and will usually take a non- 
integer value close to, but not identical with, the particle number projected or l^ ^ ^^ i 

We address here only HFB states with pairing correlations among particles of 
the same isospin. In this case, the nuclear HFB state is the direct product of separate 
HFB states for protons and neutrons, respectively, which are separately projected 
afterwards on the number of the respective particle species. 



3.1. Variational configuration mixing 

In a mean-field calculation, one often encounters a situation where the total bind- 
ing energy (mean field or projected) varies only slowly with a collective degree of 
freedom. In such a case, it can be expected that the nuclear wave function is widely 
spread around the mean-field minimum, which is beyond the scope of (projected) 
mean-field theory. These fluctuations around a single state can be incorporated 
within the generator coordinate method (GCM). The mixed projected many-body 
state is set-up as a coherent superposition of projected mean-field states \q) which 
differ in one or several collective coordinates q 

The weight function fk (q) is determined from the stationarity of the GCM ground 
state, which leads to the Hill-Wheeler-Grifhn equatiorlSl 

J2[m\q')-E,{qW)] Mq')=0, (6) 

that gives a correlated ground state, and, in addition, a spectrum of excited states 
from orthogonalization to the ground state. The weight functions fkil) are not 
orthonormal. A set of orthonormal collective wave functions in the basis of the 
states \q) is obtained from a transformation involving the square root of the norm 
kernei^. It has to be noted that projection is in fact a special case of the GCM, 
where degenerate states are mixed. The generators of the group involved define the 
collective path, and the weight functions are determined by the restored symmetry. 

For a state which would result from the mixing of different unprojected mean- 
field states, the mean particle number will not not be equal to the particle number 
of the original mean-field states anymore. Projection on particle-number, as done 
here, eliminates this problem, otherwise a constraint on the particle nu mber has to 
be added to the Hill- Wheeler- Griffin equation ([5]), see, for example, Ref.[ini. 
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The technical challenges of a configuration-mixing calculation come from the 
non-diagonal kernel s of diff erent mean-field states, which are evaluated with a gener- 
alized Wick theorerrffSn^. We represent the single-particle states on a 3-dimensional 
mesh in coordinate space using a Lagrange mesh techniqu4^. As a consequence, 
the two sets of single-particle states representing the intrinsic HFB states enter- 
ing the kerne ls are usually not equivalent, which has to be carefully taken into 
accounllSIlI!. 

We usually c omb ine the techniques presented above with a projection on an- 
gula r rnornentuirP^ that will not be discussed here. Applications are published in 
Refsl ^Q I ^^ I ^^ I ^'^ l ^^ l ^^ l ^^ i Similar methods, but without particle-number projection, 
have been set up using the Gogny force^^ and relativistic Lagrangians^S". 

For the sake of simple notation, we have introduced the GCM using a many- 
body Hamiltonian H . All the methods just mentioned, however, have in common 
that they are not based on a Hamiltonian, but an energy functional to calculate the 
binding energy. The necessary generalization will be sketched in section [5] below. 



4. Dynamical pairing correlations 

There is fundamental problem with projection after variation as performed here: it 
cannot be expected that the projection of the mean-field ground state (after vari- 
ation) gives the minimum of the energy hyper-surface that is obtained by the pro- 
jection of all possible mean-field states (which would be found by projection before 
variation). When projecting deformed mean- field states on angular- momentum, the 
intrinsic deformation of the mean-field ground state will indeed usually be different 
from the intrinsic deformation of the state giving the minimum of the projected en- 
ergy curve, with the the exception of well-deformed heavy nuclei in the rare-earth, 
actinide, and transactinide regions. 

In the context of angular-momentum projection, we overcome this problem to 
some extend by a minimization after projection (MAP), where we generate an en- 
ergy curve of projected mean-field states with different intrinsic quadrupole defor- 
mation, whose minimum provides a first order approximation to projection before 
variation. When performing a GCM calculation of projected states, the spacing of 
points along this deformation energy curve does not need to be very dense, it even 
does not have to contain the actual minimum, as the variational projected GCM 
calculation of two (non-orthogonal) projected states around a minimum has the 
ability to (implicitely) construct this minimum, as the projected state representing 
the minimum has a non-zero overlap with the actually used states. As the pro- 
jected GCM ground state also describes the energy gain from fluctuations around 
this minimum, its energy will be even below that of the minimum of the projected 
energy curve. It has to be stressed that the correlations from fluctuations around 
the projected state are outside the scope of projection before variation. As a con- 
sequence, projection before variation does not per se give a better description of 
correlations beyond the mean field than a GCM of states projected after variation. 



4, 2008 22:17 WSPC/INSTRUCTION FILE 



bender ' kazimierz 



Pairing correlations beyond the mean field 5 



While this is a standard procedure in the context of angular-momentum pro- 
jection, where a MAP is performed with respect to quadrupole deformation, it has 
been rarely addressed in connection with particle-number projection. The most ob- 
vious degree of freedom to search for a minimum of a particle-number projected 
energy curve is the amount of pairing correlations contained in the (unprojected) 
intrinsic state. Using this degree of freedom in a projected GCM calculation is 
equivalent to including the ground-state correlations from pairing vibrations (see 
Refs.Slfor an overview and Refs.^^ "^'^ for early GCM calculations using schematic 
models). Exploratory studies alon g these lines within the context of realistic mean- 



on angular momentum. 

We will present here a similar study of the role of dynamical pairing correla- 
tions in ^^"Sn. First of all, we have to note that finding a suitable constraint on 
pairing correlations is not an easy task. The obvious choice in schematic models 
is the pairing gap obtained from a pairing f orce w ith constant matrix elements. 
This coordinate was, in fact, also used in Refs.l223H_ jt has, however, some serious 
drawbacks as this pairing force leads to unrealistic asymptotics of the HFB state at 
large distances from the nucleuJ^. Constraints on other observables that measure 
pairing correlations pose similar problems even when used with a realistic pairing 
interacti on, or they introduce ambiguities on how to put them into the variational 
equations^-i ^j^q situation is similar to a constraint on a multipole moment of the 
density distribution of order i used to generate an energy surface as a function of 
deformation. Such a constraint has always to be damped at large distances from 
the nucleus, as it introduces along some direction a contribution to the constrained 
single-particle potential that diverges as —r^. 

A not completely satisfactory, but well working constraint on the amount of 
pairing correlations is provided by the strength of the pairing force. A calculation 
is done in two steps: First, the HFB or HFB-LN equations are solved using a 
"generating pairing strength" Vc- Then, in a second step, the energy of each HFB 
state is re-calculated without iterating the HFB equations using the realistic pairing 
strength Vq = —1000 MeV fm'^, either with or without projection. In a third step, 
the projected HFB states can be mixed within the GCM. Only the pairing strength 
of the neutrons is changed. For protons, we always use Vo = —1000 McV fm'^, and 
the same method as for the decription of the neutrons; in the case of pure HFB 
this means that the proton pairing breaks down. An example of such a calculation 
is shown in Fig. [1] The most remarkable findings are 

1) The HFB equations without LN corrections break down to HF at values of 
the pairing strength around —500 MeV fm'^. For those states, our formalism cannot 
introduce pairing correlations beyond the mean-field. Just above the transition from 
an unpaired to a paired system, the energy gain from projection rises very rapidly 
up to about the Vc corresponding to the minimum of the energy curve (which 
differes on the order of 10 % from Vq), and decreases more slowly afterwards. This 
is analogue to what is found in the angular-momentum projection of quadrupole 



field models are presented in Refs.l^ 




'. They do, however, not involve projection 
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Fig. 1. HFB and HFB-LN mean-field energy curves, the latter with and without the LN correction 
term in the energy, and particle- number projected energy curves as a function of the "generating 
pairing strength" Vc for a spherical mean-field state in ^^"Sn. The dots with horizontal bars 
represent the energy of the projected GCM ground state plotted at its average "generating pairing 
strength" . The energy scale on the left is normalized with respect to the HF ground state, while 
the scale to the right gives the total binding energy. 



deformed states around a spherical configuration. 

2) There are obvious differences between HFB and HFB-LN energy curves, 
both on the mean-field level and projected, for small generating pairing strength 
Vc < —1500 MeV fm^, which is indeed in the regime of the physical pairing strength. 
In the strong pairing regime at larger values of Vc, the energy curves obtained from 
projection of the HFB and HFB-LN states are nearly identical. The remaining 
difference might be attributed to the presence of proton pairing correlations in the 
HFB-LN case, while they are absent for HFB. 

3) The LN correction overestimates the energy gain from projection at small 
generating pairing strength, and underestimates it in the strong pairing regime. 

4) The additional energy gain from the GCM of projected states is quite small, 
around 100 keV when starting with HFB-LN states, and about 200 keV for pure 
HFB states. The unphysical breakdown of HFB in the weak pairing regime leads 
to a smaller binding energy at the minimum of the projected energy curve, and 
gives a potential energy surface that is stiffer than the HFB-LN one for small Vc- 
This leads to a GCM ground state from projected HFB that is less bound than 
the GCM gound state from projected HFB-LN. This also pushes the GCM ground 
state wave function from projected HFB to larger values of Vc than the one from 
projected HFB-LN, as seen from the larger average generating pairing strength 

= 'l2vc 9oi^c)Vc in the projected HFB case. 

This calculation, of course, scratches only on the surface of the importance 

of dynamical pairing correlations. The question of better constraints on pairing 
correlations, systematics, excited states, and the coupling with deformation modes 
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Fig. 2. Lower panel, left: particle-number projected binding energy of calculated with the 
Skyrme interaction Sly4 and a density-dependent pairing interaction as a function of the mass 
quadrupole deformation for 9 and 99 discretization points for the integral over the gauge angle (/3 
in Eq. l(4]l. Upper panel, left: canonical single-particle energies (full lines: parity dotted lines: 
parity —1) and Fermi energy (dashed line) for neutrons. Right panel: dispersion of the neutron 
number (A^o l-'V^ | A^o) — (-'Vq |7V|7Vo)^ for 1 (no projection), 3 and 5 discretization points. 



will be addressed elsewhere. 



5. The divergence in particle-number projection 
5.1. General features 

It has been noticed for a long time that the particle-number projected energy might 
exhibit divergenceJ23ES] -^^^ej^ a single-particle level crosses the Fermi energy and 
has the occupation = 0.5. More recently, it was pointed out by Stoitsov et aM^, 
that in addition to the divergence there appears a finite step in the projected energy 
when passing this situation as a function of a constraint. An example is given in 
Fig. [21 where two clear divergences appear when a fine discretization of the integral 
over gauge angles, Eq. (g]), is used. The steps are also hinted, with the binding 
energy changing on both sides from a lower curve around sphericity to a higher 
lying curve when passing the divergence. 

First of all, it has to be stressed that the binding energy is the only observable 
that shows an anomaly when a single-particle level crosses the Fermi energy. Neither 
the overlap, nor any observable calculated from any n-body operator shows an 
unusual behaviour. This is examplified in Fig. [5] by the dispersion of the particle 
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number, a two-body operator that also provides a measure for the numerical quality 
of projection. For a small system like ^^O, the integrals over gauge angles are 
already numerically converged for L = 5 discretization points. The only exception 
is the binding energy. However, an extremely huge number of discretization points 
is needed to see the divergence develop in (and even more for heavier nuclei), 
which is one of the reasons why the divergence often remains undetected. The other 
reason is that it is very unlikely that one of the discrete points used to calculate a 
potential energy surface hits the narrow region where the divergence appears. This 
is different when performing projection before variation, where the variation will 
detect the divergences more easilj^^Sl. 

The appearance of a divergence for the binding energy, but no other observable 
is related to the the particular definition of the binding energy in our method: the 
energy is calculated from an energy density functional, while everything else is, for 
the moment at least, calculated as the expectation value of an operator. 

To understand the origin of the divergence, we have to look into the defini- 
tion and evaluation of the energy functional. For the sake of simple notation and 
a transparent argument, we will restrict ourselves here to a toy model with one 
kind of particles only, and a two-body interaction. The additional complications 
introduced by dens ity dependencies as needed in realistic energy functionals will be 
discussed elsewhere^. The further generalization to a system composed of protons 
and neutrons is then straightforward. All expressions given below are evaluated in 
the canonical basis, as this basis will turn out to be the only basis in which the 
origin of the divergence can be clearly identified and analyzed. 



5.2. The Hamiltonian case 

As a reference, for which everything is properly defined and no problem occurs, we 
will use the energy obtained from a two-body force. At the HFB level, one has 

£"[p, = fkinb] + i ^'J™" (HFB|ata]a„a„|HFB) (7) 

The contractions are the usual density matrix and pair tensor in the canonical basis 
/ t \ (HFB|ata,jHFB) c 2 c / s 



^ - \ 01 (HFB|a.a^|HFB) e . /qs 

a.flp) = 't^. = (HFB|HFB) " " "^"^ 



10 (HFB|«t«t|HFB) ^ _ 

Unlike in papers on standard HFB theory we distinguish already here between two 
different pair tensors, as they generalize differently for particle-number projected 
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HFB states. Unless necessary, we will not specify the kinetic energy in what follows. 
It is given by a one-body operator, always evaluated as such from whatever the left 
and right many-body states might be, and free of any divergence problems. The 
generalization of £^ to particle-number projected states is straightforward 

^ (HFB|gP^|HFB) ^ f'''^^ (HFB(0)|i/|HFB(^)) (11) 

(HFB|P^|HFB) Jo SvrPjVo 

with Vpfg — (HFB|P^|HFB). The relevant piece is the calculation of the Hamilto- 
nian kernel (HFBj/f |HFB((/5)) Although for particle-number projection the Hamil- 
tonian kernel can - in principle - be evaluated using the stand ard Wi ck theorem, 
it is much more convenient to apply a generalized Wick theorenPHZl 



(HFB 1 7^1 HFB (^)) = 



Av)- (12) 



The basic contractions are given by 



y2 g2i¥> 

(at a^)y = PtJ.v{^p) = '^i^M = 2 r 2 ^2»y ^''t^ (12) 

U,, -f- v.. 6 



UnVr, 



1 , 1 ^ 



and the norm kernel T{(p) by 



X(^) = (HFB|HFB(^)) = n + < ^'"')- (16) 



5.3. The EDF case 

In the case of general energy density functional (EDF) in the spirit of a Kohn-Sham 
approach with pairing, the energy is given by 

= £ki„[p] +£pp[p] (17) 

where fkin is the kinetic energy, Spp the energy of the particle-hole interaction, and 
£kk the energy from the particle-particle (pairing) interaction. The only restriction 
that we impose on the energy functional for the moment is that it is bilinear in 
either the density matrix or the pair tensor, which is done to keep the analogy with 
the Hamiltonian case. Then, the energy functional can always be described in terms 
of the kinetic energy and a double sum over not antisymmetrized two-body matrix 
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elements wPp^j^ and w"", different in the particle- hole 

= £kin[p] + Y («|.«/^) + E <A>^P(«|.4)(«'^«'')- (18) 



This is the Kohn-Sham approach to an EDP^, formally generalized to systems with 
pairing by Oliveira et al^^^ where the actual density (matrix) is calculated from an 
auxiliary independent quasi-particle state |HFB). Equation (|18p might appear to be 
an unusual representation of an energy functional, but, for example, the interaction 
of a term bilinear in the local density, where f{\r — r'\) is an arbitrary function and 
C the coupling constant, translates as 



C jjdPr d^r' p{r) /(|r - r') p{r') 

= ^ E [[ d'r d^r' p,,i;lir)4,^ir) fi\r ~ r'\) p,,^lir')Mr') 



= C ^ (19) 

Similar expressions are obtained for any other bilinear contribution to an energy 
functional that does not have density-dependent coupling constant. 

For the generalization of the energy functional to the case of particle-number 
projected states, the same generalized Wick theorem as above is used to define 

f™F(^) - r,^^ e-^^^" H™F(^) (20) 



where the Hamiltonian kernels are now given by 

(21) 

The identification of the origin of the divergence proceeds as follows 

1) As analyzed in detail by Anguiano et al\^, the divergence appears for those 
terms in the energy that originate from the interaction of a particle with its conju- 
gated partner 



H^"^(^) = 



„,2 2iip ,,2p2iifi 
PP , PP , PP \ P P 

111- - - - + - -\- ?v; - - 1 — 

'^PPPP ' "^Pf^PP ^ '^ppppl y2 _|_ ^2g2ii^ y2 ^_ y2^2if> 



' ^'-^ II II II II 



I{ip). (22) 



For — 0.5 and ip = 7r/2, the denominator in the transition densities becomes 
zero. One of the two denominators is canceled by the same factor contained in 
I{(p), Eq. Uni), while the other one causes the divergencJ 
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2) It can be shown from very general argument that the Hamiltonian kernel 
should have a certain dependence on (p. The divergent contributions to (l22|) do not 
follow this rule. This shows that such terms are spurious even for ip ^ 7r/2. 

3) In the case of a two-body Hamiltonian the divergence disappears when one 
identifies 2w^p^^ = iw'^p^p^ = w^ip^p and w^^^^ — 0. This can be used to combine 



the and w^e^"^ factors such that they cancel the dangerous denominator 

4) The matrix elements of the kind it'^^^^ might indeed have non-zero values 
in an EDF, which represent a spurious interaction of a particle with itself. They 
violate the exchange symmetry in a Fermionic systems (the Pauli principle), and 
lead to a spurious contribution to the total binding energy already on the mean-field 
level. The appearance of this so-called " self-interaction" is a well-known annoyance 



of EDFs for electronic systema^Ql, but ignored in nuclear physics so far. 

5) In the case of an EDF with pairing, there is an additional spurious "self- 
pairing" interaction that originates from the scattering of a pair onto itself. The 
interaction energy from two isolated Fermions occupying pair-conjugated orbitals 
divided by the occupation of the pair should have the same value as if pairing 
correlations were not presenile. Again, this is violated by a general EDF. The 
actual expression for the spurious self-pairing energy combines matrix elements 
and the occupation factors that weight them and will be given elsewher^^. 

6) The divergent contributions to Eq. ([22| come from terms that represent self- 
interaction and self-pairing. By that, the projection adds a second level of spuriosity 
to these terms. Not the whole contribution from the self-interaction and self-pairing 
terms is divergent, though. 

7) The self-interaction and self-pairing contribute to the total energy not only 
because the relevant matrix elements have spurious non-zero values, but also be- 
cause they are multiplied with unphysical weights. When evaluating the Hamilto- 
nian kernel (|12p by commutating the creation and annihilation operators of the 
Hamiltonian with those setting up the HFB states until they hit the vacuum, it 
becomes clear that each pair of conjugated particles (/x, /x) can be multiplied with 
one e^"'' factor stemming from the rotated state to the right only, and that the 
possible combinations of and for a given fi will always be bilinear. Only when 
the Wick theorem is applied, one obtains terms which are of 4th order in and 

and quadratic in e'^'''^ . The reason is that multiple contributions from the same 
particle or pair of particles are not excluded from the sum as the Wick theorem im- 
plicitely assumes that the matrix elements these terms multiply are zero or sum up 
to zero, so that they will not contribute anyway. In an EDF these matrix elements 
are not zero anymore, but the Wick theorem is still used. 

One rigorous way to remove the spurious terms from an energy functional would 
be to remove all possible self-interactions in the total energy by excluding that 
the same summation index appears more than once in a contribution to Eq. (jl8p 
This w as, in fact, already used by Hartree in his seminal paper on the Hartree 
method^. Or, alternatively, one sums up explicitely all the contributions where 
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the same summation index appears two times or more often and subtracts this as a 
self-energy correction from the total energj^^. This leads, however, to enourmous 
complications in the variational equations, particularly when extended to density- 
dependent interactions. 

5.4. The origin of the finite step 

As pointed out by Stoitsov et a/I^Sl, there is also a finite step that appears when 
passing the divergence as a function of a collective coordinate. The analysis of its 
origin will be given elsewherj22l. Let us just outline the main arguments: With the 
substitution z = e*"^, the integral over the real gauge angle ip in Eq. ([2]) can be 
transformed into an integ ral in the complex plane, that can be analyzed with the 
tools from function theorjl2Si32| jj^ ^^le Hamiltonian case, the integral has a pole at 
z = 0, which has an order of the number of particles below the Fermi energy. The 
residue of this pole is proportional to the projected energy. In a projected EDF 
framework, there are two changes to this scenario: First, there is an additional 
contribution to the pole at z = 0, this time from all particles, and second, there 
appear additional poles at = ±i|m^|/|i;^| along the imaginary axis, which also 
contribute for all particles below the Fermi energy. Both originate from the same 
terms as the divergence. The two additional contributions to the projected energy 
are huge (on the order of several hundered MeV for ^^O), but of opposite sign 
and nearly, but not exactly, canceling each other. The step appears when one of 
the poles at ±i|m^|/|z;^| enters or leaves the integration contour |z| = 1 when the 
single-particle energy of the corresponding particle crosses the Fermi energy. 

6. Further discussion 

Density-dependent terms add further complications to the divergence, which we will 
not address here in detail. Let us just comment that, first, in a projected theory, 
the evaluation of a density-dependent term with non-integer power requires the 
evaluation of a (multivalued) root of a complex number, which leaves the Riemann 
sheet for choice. Second, when one expands the densities in the energy functional for 
density-dependent terms, and combines the resulting terms similar to Eq. (|2ip . one 
finds again terms which contain more than just one v^e^^^ factor originating from 
|IIFB((^)) (one of them with a usually non-integer power), which again introduces 
an unphysical dependence on (f into the Hamiltonian kernel. 

There are good and profound reasons to use different effective interactions in the 
particle-hole and particle-particle channels. The particle-hole and particle-particle 
channel of the effective interaction sum up different classes of diagrams^, which 
inevitably leads to different expressions for the effective interaction in both chan- 
nels. Short-range correlations are resummcd into the functional providing the two 
channels with different density-dependent forms. On the other hand, the long-range 
correlations from large-amplitude fluctuations around the mean-field states are ei- 
ther neglected in a pure mean-field approach when they are assumed to be small, or 
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otherwise described explicitely by projection and GCM-type configuration mixing. 
A correction that removes the divergent part and the finite step from the projected 
energy can be set-up when identifying the terms with an unphysical dependence 
on the gauge angle com paring the expressions obtained from the standard and 



While particle-number projection is the prominent example for the appearance 
of a divergence, similar divergences can be expected in any GCM calculation or 
projection on any other quantum number. As in these cases the mixed states have 
different canonical bases, the transition density matrix is usually not diagonal and 
the analysis is less obvious. The widely used collective Schrodinger equation and 
Bohr Hamiltonian that are set-up through a series of approximations (at some 
points including improvements) from exact configuration mixing cannot be expected 
to be free of problems stemming from the divergence either, although there they 
will be more difficult to identify. 

7. Summary 

First, we have examined for the example of ^^"^Sn the effect of a minimization after 
particle-number projection and of GCM ground state correlations from pairing 
vibrations. For ^^"Sn, the effect is visible, but not enormous, and clearly smaller 
than the current uncertainties from the parameterization of the effective pairing 
interaction. This might be, however, different for other nuclei. 

Second, we have examined the divergence that appears in particle-number pro- 
jected EDF. As its origin we identify the spurious self-pairing contribution con- 
tained in all common energy functionals for self-consistent mean-field models. On 
the mean-field level, the self-pairing as such adds a small spurious contribution to 
the total binding energy, that is similar to, but smaller, than the usual spurious 
centre-of-mass or rotational energies from broken symmetries (in fact, together with 
the self-energy it is the spurious energy from violating the Pauli principle in approx- 
imate EDF approaches). When generalizing the EDF to particle number projection 
via the generalized Wick theorem, the self-pairing provides the Hamiltonian ker- 
nels with an unphysical dependence on the gauge angle, which ultimately leads to 
divergences and steps. 
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